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Chapter  I 


INTRODUCTION 


1.1  INTRODUCTION 

A  free-edge  in  a  composite  is  the  intersection  of  an  interface 
plane  (between  any  two  layers)  and  the  free  surface  of  the  composite, 
while  a  contact-edge  is  the  intersection  of  an  interface  plane  (between 
any  two  layers)  and  the  contact  surface  of  the  composite.  The  unusually 
large  and  possibly  infinte  stress  at  the  free-edge  and  the  contact-edge 
is  one  of  the  factors  responsible  for  delamination  when  the  composite  is 
subjected  to  external  loading.  Since  more  composites  are  now  used  in 
space  vehicles,  it  is  important  to  analyze  the  nature  of  the  stress 
singularities  at  the  free-edge  and  the  contact-edge  so  that  special 
attention  can  be  paid  to  the  design  of  mechanical  fastenings  and  joints. 
The  research  reported  here  is  part  of  the  effort  of  the  "predictive 
response  capability  development"  program. 

There  have  been  many  investigators  who  analyzed  the  stress  near  the 
free-edge  [1-9].  An  analytical  solution  which  is  valid  for  the  whole 
composite  is  practically  impossible  to  obtain.  Several  approximate 
numerical  solutions  are  available  which  show  a  good  agreement  between 
each  other  for  points  away  from  the  free-edge.  For  points  near  the 
free-edge,  numerical  solutions  are  not  capable  of  predicting  an  infinite 
stress  when  it  exists,  and  this  is  where  the  discrepancies  between 
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various  approximate  solutions  occur.  Wang  and  Choi  [8]  used  an 
eigenfunction  expansion  technique  to  determine  the  stress  in  the 
interface.  However,  the  completeness  of  the  eigenfunction  expansion  is 
an  open  question  [10].  In  fact,  the  existence  of  the  logarithmic 
singularity  discussed  in  this  paper  implies  that  the  eigenfunction 

expansion  in  terms  of  r^  powers  may  not  be  complete.  It  is  doubtful 
that  the  addition  of  (In  r)  terms  would  make  the  eigenfunction  expansion 

2 

complete.  As  pointed  out  in  [10,  11],  singular  terms  of  (In  r)  and 

3 

(In  r)  etc.,  may  also  exist. 

For  composites  whose  layers  are  isotropic  elastic  materials,  use  of 
the  biharmonic  function,  or  the  Airy  stress  function,  seems  to  be  the 
universal  approach  in  the  analysis  of  the  stress  singularities.  (See 
[10,  12,  13],  for  example).  There  appears  to  be  no  universal  approach 
in  analysing  the  stress  singularities  in  anisotropic  elastic  materials. 
Lekhnitskii  [14]  introduced  two  stress  functions  to  analyze  general 
anisotropic  materials.  His  approach  was  used  by  Wang  and  Choi  [8]  to 
study  the  thermal  stresses  at  the  interface  in  a  layered  composite. 

Green  and  Zerna  [15]  employed  a  complex  function  representation  for  the 
general  solution.  Their  approach  was  used  by  Bogy  [16]  and  Kuo  and  Bogy 
[17]  in  conjunction  with  a  generalized  Mellin  transform  to  analyze 
stress  singularities  in  an  anisotropic  wedge.  In  this  paper  we  use  the 
approach  which  was  originated  by  Stroh  [18]  and  further  developed  by 
Barnett  and  others  [19*21]  for  studying  the  surface  waves  in  anisotropic 


elastic  materials. 
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*  5 

While  the  nature  of  the  singularity,  be  it  k  r  or  k(ln  r) ,  is 
independent  of  the  stacking  sequence  of  the  layers  in  the  composite  and 

JU 

the  complete  boundary  conditions,  the  unknown  constant  k  in  the 
singular  solution  is  not.  This  suggests  that  one  might  use  a  special 
finite  element  at  the  free-edge  or  at  the  contact-edge  (with  regular 
finite  elements  elsewhere)  so  that  the  exact  nature  of  the  singularity 
is  prescribed  in  the  special  element  while  the  unknown  constants 
associated  with  this  special  element  are  determined  by  solving  the 

complete  boundary  value  problem.  If  k  so  obtained  happens  to  be  zero 

at  a  particular  free-edge  or  contact-edge,  there  is  no  r^  singularity 
at  that  point.  The  purpose  of  this  paper  is  to  provide  the  exact  nature 
of  the  singularity  at  the  free-edge  and  at  the  contact-edge. 

In  the  following  sections  of  Chapter  I,  the  notation  for  the  basic 
equations  (strain-displacement,  stress-strain  and  equilibrium  equations) 
is  explained,  and  expressions  for  the  elasticity  constants  and  elastic 

•S  A 

compliances  of  materials  orthotropic  in  the  (x^^,  x^,  x^)  coordinate 
system  are  given. 

In  Chapter  II  a  free-edge  problem  which  undergoes  a  two-dimensional 
plane  strain  deformation  is  analyzed.  The  assumed  expressions  for 
stress  and  displacement  when  applied  to  the  boundary  and  continuity 
conditions  result  in  a  system  of  linear  homogeneous  equations.  The 
eigenvalues  6  of  this  system,  both  real  and  complex  are  listed,  and  the 
procedure  used  to  find  the  complex  values  of  6  is  explained.  The 

existence  of  a  singularity  of  k  r^  for  negative  6  is  discussed. 


u 


In  Chapters  III  and  IV  the  free-edge  problem  is  analyzed  for  three- 
dimensional  deformation.  In  Chapter  III  it  is  established  that  the 
superposition  of  the  uniform  extension  term  on  the  two-dimensional 

problem  of  Chapter  II  is  a  first  order  approximation  to  a  three- 
dimensional  deformation.  In  Chapter  IV  it  is  shown  that  the  addition  of 
the  uniform  extension  term  to  the  equations  of  Chapter  II  results  in  a 
system  of  non-homogeneous  linear  equations.  Two  methods  are  then  used 
to  solve  for  the  assumed  homogeneous  stress.  These  solutions  are  valid 
for  (0/-0)  and  (0/90)  composites.  A  modified  method  is  then  introduced 
to  solve  for  the  non-homogeneous  stress  of  the  other  (0/0') 
combinations.  The  existence  of  the  .k(ln  r)  singularity  is  discussed, 
and  the  logarithmic  stress  intensity  factor  k  is  listed  for  various 
combinations  of  (0/0').  Values  of  which  depend  on  the  angle  0  in 

polar  coordinates  are  also  listed. 

In  Chapter  V  a  contact  problem  is  analyzed.  The  altered  boundary 
conditions  lead  to  a  new  system  of  linear  homogeneous  equations  for  the 
two-dimensional  plane  strain  problem.  The  roots  6,  both  real  and 

complex,  are  listed  and  the  existence  of  the  singularity  k  r°  is 
discussed.  The  three-dimensional  deformation  problem  is  then  analyzed, 
and  a  particular  solution  which  is  uniform  in  stress  is  found  and  listed 
for  some  (0/0')  combinations. 

1.2  BASIC  EQUATIONS 

In  a  fixed  rectangular  coordinate  system  x^,  (i  =  1,  2,  3),  let 
Uf’  ®ij  ®ij  displacement,  stress,  and  strain,  respectively. 
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The  strain-displacement,  stress-strain  and  equilibrium  equations  may  be 
written  as 


E  .  .  =  (3u./3x.  +  3u./3x. )/2 
ij  1  J  J  1 


°ij  ^'ijkt^kH 


(1.1) 

(1.2) 


£  .  .  =  S  .  -O,  , 

ij  ijki  kH 


30  .  ./3x.  =  0 
ij  J 


where  repeated  indices  imply  summation,  and 


^ijke  ~  ‘^kiij  ‘^jike 


®ijk£  ~  ^kilij  ^jikfc 


(1.3) 

(1.4) 

(1.5) 

(1.6) 


are  the  elasticity  constants  and  the  elastic  compliances,  respectively. 
Due  to  the  symmetric  property  of  Eqs.  (1.5,  1.6),  one  may  rewrite 


(1.7) 

(1.8) 


where 


1.5, 

,  1.6) 

as 

0  = 

c  , 

c 

=  c 

q 

qt  t 

qt 

tq 

e  = 

s 

=  s 

q 

qt  t* 

qt 

tq 

0 .  = 

0 

0  = 

0  = 

o„_ 

1 

11’ 

2 

22’ 

3 

33 

Q 

II 

®23’ 

°5  = 

‘’13’ 

‘’e  = 

='12 

=  1  = 

=  11’ 

=  2  = 

=  22’ 

=  3  = 

=  33 

^  = 

^=23’ 

=  5  = 

^=13’ 

=  6  = 

^=12 

(1.9) 


(1.10) 
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The  transformation  between  c,.,  .  and  c  are  as  follows 

ijki  qt 


ij  or  kH  = 


q  or  t  = 


(1.11) 


Since  the  transformation  between  a.,  and  a.  is  not  identical 

ij  1 

to  that  between  e . .  and  e . ,  in  order  to  represent  the  transformation 
1  j  1 

between  the  subscripts  of  s,.,  ,  and  s  ,  Eq.  (1.11)  must  be  modified 

ijki  qt 

as  follows.  If  either  q  or  t  is  larger  than  3,  s  ^  =  2s..,..  If 

^  **  qt  ijkH 

both  q  and  t  are  larger  than  3,  s^^  =  4s  [22]  contains  further 

discussion  about  these  transformations. 


1.3  ORTHOTROPIC  MATERIALS 

For  orthotropic  materials  with  the  (x^^,  x^ .  x^)  axes  of 
symmetry,  s^^  are  zero  except  [25,  26) 

Sii  =  1/E^,  S22  =  I/E2,  5^3  =  l/E^,  1 


®44  ~ 

I/G23, 

"55  "  ^^Si’ 

"66  =  1/^12 

"12  = 

"21  = 

V21/E2 

^  (1.12) 

"13  = 

"31  = 

V31/E2 

®23  “ 

S32  = 

V32/E3 

4 
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where  E^.  E^,  E3.  G23.  63^,  '^21*  '^31*  ^32 

the  engineering  constants.  The  relation  between  s^j  and  may  be 

found  in  [22,  25,  26], 

Consider  a  laminated  composite  which  consists  of  a  finite  number  of 
anisotropic  elastic  layers  perfectly  bonded  at  the  interface.  Fig.  1. 
Each  layer  of  the  composite  lies  in  a  plane  parallel  to  the  (x^^,  X3) 

plane.  A  cross  section  of  the  composite  is  shown  in  Fig.  2.  This 
angle-ply  graphite/epoxy  laminated  composite  is  assumed  to  be 
orthotropic  with  respect  to  the  (Xj^,  X2,  X3)  coordinate  system. 

The  X3-axis,  which  is  the  direction  of  the  fibers  (as  shown  in  Fig. 


3),  makes  an  angle  0  with  the  X3-axis.  It  will  be  assumed  that 

materials  1  and  2  in  Fig.  2  are  made  of  the  same  orthotropic  material, 

but  that  the  orientation  of  the  fibers  is  different;  the  angles  9  and  0' 

vary  between  layers.  Therefore,  c..  and  s..  are  the  same  for  both 

ij  ij 

materials,  while  c..,  s,.  may  not  be.  The  relation  between  c,, 
ij  ij  ij 

and  c. .  may  be  found  in  [221.  The  relations  between  s..  and  s.. 
iJ  ij  ij 

may  be  obtained  from  the  relations  between  c^^  and  c^^  by  replacing 
c^j  by  s^j  with  the  following  exceptions.  If  either  i  or  j  is 


larger  than  3,  replace  2c^j  by  If  both  i  and  j  are  larger  than 

3,  replace  by  The  same  rules  apply  in  replacing  c^^  by 


s  .  .  . 
^  J 


Therefore  we  have 


^11  =  Sjj  +  c  s  (2Sj3  +  S33)  +  s  S33 


2-  2- 
^12  =  '=  "12  ^  "  ^32 


4  .  4 


2  2/^  ,  ^  c: 


~  (c  +  s  )s^3  +  c  s  ^33  ^55^ 


Si5  =  cs[c2(2s^3  +  Sgg  -  2ajp  -  s  (2^3^  +  S33  -  2S33)] 


®22  ~  ®22 


2-^  2/n 

®23  =  ‘^  ^23  "  ^21 


*25  ~  2cs(s23  -  s^j) 


(1.13) 


Aa  22,/ijA  /*V  ..  4y\ 

*33  =  *33  *  ^2=13  +  S33)  +  s  Sj 


S35  =  '^s[-c  (2^3^  +  S55  -  2S33)  +  s  (2s^3  +  §33  -  2S^^)] 

2^  ^  2^ 

®44  ~  ^  *44  *  *66 

*46  ■  ^^*^*44  ■  *66^ 


S55  -  (c  -  s  )  S33  +  4c  s  (s^j  -  2s^3  +  S33) 


2^  ,  2^ 
s,,  =  c  s,,  +  s  s.. 
66  66  44 


where,  for  simplicity,  the  notation  c  =  cos9,  s  -  sinB  has  been  used. 


All  other  terms  of  s^^  are  zero. 


Chapter  II 

TWO  DIMENSIONAL  FREE-EDGE  SINGULARITY  ANALYSIS 


2.1  FORMULATION  OF  TTO  MATRIX  K 

In  order  to  study  the  stress  at  a  free-edge  of  an  interface  of  a 
laminated  composite,  the  origin  of  the  (x^^,  X2,  x^)  axes  is  placed 

at  one  of  the  free-edge  points,  as  shown  in  Fig.  2.  The  (x^^,  x^) 

plane  is  along  the  interface,  while  the  (X2,  x^)  plane  is  the  free- 
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where 


°ik  =  Qik  ^  P^^ik  ^  ^  P  ^ik 


(2.4c) 


^ik  ‘^ilkl’  ^ik  ~  ‘^ilkZ’  ^ik  ~  ‘^i2k2 


(2.4d) 


For  a  non-trivial  solution  of  u^,  it  follows  from  Eq.  (2.4a)  that  the 

determinant  of  must  vanish.  This  results  in  a  sextic  equation  for 

p.  Since  the  eigenvalues  p  are  all  non-real  [22,  14,  18],  there  are 
three  pairs  of  complex  conjugates  for  p  and  three  pairs  of  associated 
eigenvectors  u^.  For  isotropic  materials,  all  p's  have  the  value  +i  . 

To  analyze  the  singular  nature  of  the  stresses  at  the  origin,  the 
function  f(Z)  is  chosen  to  be 


f(Z)  =  Z^  ®/(l  +  6) 


(2.5) 


where  6  is  a  constant.  Equations  (2.1,  2.2)  for  displacement  and  stress 
can  then  be  written  as 


(2.6) 


0..  =  Z  (A-T..  -Z/  +  B-T..  ,z/} 
1 j  L  ij ,L  L  L  ij ,L  L 


(2.7) 


where  an  overbar  denotes  the  complex  conjugate;  and  B,.  are 

Ij  Li 

complex  constants,  and  the  subscript  L  identifies  the  three  pairs  of 
eigenvalues.  Unless  otherwise  indicated,  Z  in  Eqs .  (2.6,  2.7)  and  in 
the  sequel,  stands  for  summation  over  L  from  L  =  1  to  3.  Using  the 
polar  coordinates  (r,  0),  Fig.  2,  Z  may  be  rewritten  as  [11,  12] 


Z  =  +  px^  =  rC 


(2.8) 


where 


i  =  cos^  +  psin^ 

Equations  (2.6,  2.7)  can  then  be  written  as 


(2.9) 
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°ij  "  ^  ^Vij.L^L  Vij.L^L  ^ 


(2.10) 

(2.11) 


Similar  equations  may  be  written  for  the  material  with  elasticity 
constants  '  by  adding  a  prime  to  all  quantities  except  r,  0  and  6. 

We  see  from  Eq.  (2.11)  that  if  the  real  part  of  S  is  negative,  is 

singular  at  r  =  0.  By  applying  the  stress  free  boundary  conditions 


Oi  =  05  =  06  =  0 


(2.12) 


at  0  =  ±ir/2  and  the  interface  continuity  conditions  at  ^  =  0 

[Uj]  =  [u^]  =  [u^l  =  0  (2.13a) 

IO2I  =  [O4]  =  l%]  =  0  (2.13b) 

where  (f)  =  f  -  f'  represents  the  difference  in  f  values  across  the 
interface,  one  obtains  12  linear  homogeneous  equations  for  Aj^, 

^L*  ’  ®l'  written  as 

K  (6)q  =  0  (2.14) 

where  is  a  complex  valued  square  matrix  whose  elements  depend  on  6, 
and  q  is  a  column  matrix  whose  elements  are  Aj^,  Bj^,  Aj^'  ,  B^^'  , 

(L  =  1,  2,  3).  If  6,  which  can  be  real  of  complex,  is  a  root  of  the 
determinant 

l|K^(6)ll=0  (2.15) 

then  a  nontrivial  solution  exists  for  q,  and  hence  for  the  stress  and 
displacement . 

When  the  root  6  of  Eq.  (2.15)  is  real  we  may  choose 
®L  =  \  =  ^ 


(2.16) 
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where  and  are  real.  Equations  (2.10)  and  (2.11)  then  have  the 
real  expressions 

°ij  =  ^  ^  V'"("ij.L'^L^)>  (2.18) 

where  Re  and  Im  stand  for  real  and  imaginary,  respectively.  Equation 
(2.14)  is  then  replaced  by 

K(6)a  =  0  (2.19) 

where  g  is  now  a  real  valued  square  matrix  and  a  is  a  real  column  matrix 

whose  elements  are  a  ,  a^ ,  a  '  and  a^ ' ,  (L  =  1,  2,  3).  The  real 

Li  1j  Li  Ij 

root  6  is  then  obtained  from 

llg(5)!l  =  0  (2.20) 

Before  we  present  a  numerical  procedure  for  finding  a  complex  root 
from  Eq.  (2.15),  we  point  out  that  6  =  0  is  always  a  root  of  Eq.  (2.15) 
and  Eq.  (2.20).  It  should  also  be  pointed  out  that  the  formulation  here 
tacitly  assumed  that  the  eigenvalues  p  of  the  elasticity  constants, 

Jj 

Eq.  (2.4a)  are  distinct.  For  degenerate  cases  in  which  p^  has  a 

Li 

multiple  root,  Eqs .  (2.6,  2.7)  have  different  expressions.  The  correct 
expressions  for  Eqs.  (2.6,  2.7)  when  Pj^  has  a  multiple  root  were  given 

in  [11]  when  f(Z)  assumes  the  special  form  of  Eq.  (2.5)  and  in  [23]  when 
f(Z)  is  arbitrary. 

2.2  PROCEDURE  TO  FIND  THE  ROOTS  OF  HK^H 

To  find  the  roots  6  of  Eq.  (2.15),  which  may  be  complex,  the 
following  method  was  used.  Let 

S  =  5  +  i»l  (2.21) 


_ 
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llK^II  =  u  +  iv  (2.22) 

The  complex  plane  5  is  divided  into  square  regions  of  prescribed  size. 
The  nodal  points  of  a  typical  square  region  have  the  coordinates 
Cjj.  -  1>  2,  3,  4),  as  shown  in  Fig.  4.  The  value  of  IlK^H  at 

each  node  of  the  square  is 

=  I-CIS.IDIj  =  5^  (2.24) 

If  the  square  contains  the  curve  u  =  0,  then  either  the  value  at  one  of 
the  nodes  will  be  zero  (in  this  case  the  curve  intersects  the  node),  or 
the  sign  of  u  will  change  between  at  least  one  pair  of  adjacent  nodes, 
Fig  5.  A  possible  exception  to  this  is  the  case  shown  in  Fig  6,  which 
may  be  overcome  by  initially  selecting  a  small  enough  square.  To  test 
if  the  line  u  =  0  is  present,  the  computer  program  checks  the  nodal 
values  Uj^  for  zero,  and  also  checks  for  a  change  of  sign.  If  either 

of  these  conditions  is  satisfied,  the  program  repeats  the  same  procedure 
to  determine  if  the  v  =  0  curve  exists  in  the  square.  When  the  program 
determines  that  both  curves  are  contained  in  the  square,  then  the  square 
is  divided  into  four  smaller  squares,  each  of  which  is  analyzed  in  the 
same  way,  and  again  subdivided  if  necessary.  Fig.  5.  A  square  which  is 
found  to  contain  only  one  of  the  curves,  or  neither  of  the  curves  is 
discarded.  The  procedure  continues  until  either  no  squares  remain,  or 
until  the  dimension  of  the  square  is  smaller  than  a  specified  error 
parameter.  If  this  parameter  is  small  enough,  then  the  fact  that  the 
two  curves  exist  in  the  square,  would  suggest  that  they  intersect,  and 
that  a  root  exists  inside  that  square.  A  flow  chart  of  the  subroutine 
can  be  found  in  Fig.  7. 
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It  should  be  pointed  out  that  for  a  very  small  square  which 
contains  the  root  6,  1|K^11  can  be  quite  large  at  the  four  nodes  of  the 

square.  For  example,  Table  1  lists  the  values  of  6  and  |lK  ||  found  by 


the  above  method  with  the  dimension  of  the  square  being  reduced  to  10 
The  values  of  at  the  four  nodes  are  in  the  order  of  10^. 


-9 


2.3  THE  ROOTS  6 

Two  different  laminated  composites  are  used  for  the  numerical 
calculations.  Each  layer  of  a  composite  is  assumed  to  be  made  of  the 
same  orthotropic  material.  The  orientation  of  the  axes  of  symmetry 

A  A 

(Xj^,  x^,  x^),  however,  differs  from  layer  to  layer.  The  following 

engineering  constants  for  the  layers  in  the  two  composites  are  taken 
from  [6,  24],  respectively. 


Composite  W 

(Typical  high  modulus  graphite-epoxy,  [6]) 


^2  = 

2.1  X  10 

psi 

s  = 

20  X 

10^  psi 

CM 

"  ^3 

=  Si  =  ■ 

85  X  10^  psi 

''21 

II 

< 

h-* 

■  ''32  ”  ■ 

21 

(2.25) 


.21 
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Composite  T 

(T300/5208  graphite  epoxy,  [24]) 


=  1.54  X  10^  psi 


E^  =  22  X  10  psi 


°12  “  °23  ~  Si  ~ 


'^21  '^31  ''32  ~ 


(2.26) 


Using  Eq.  (2.25)  or  Eq.  (2,26),  s^^  is  obtained  from  Eq.  (1.12), 

A  ^ 

while  is  computed  from  s^^  using  the  relations  derived  in  [25, 

26]  .  s  .  .  ,  s  .  . ' ,  c  .  .  and  c . . '  associated  with  various  9  and  0 ’ 
ij  ij  ij  ij 

are  then  determined,  from  Eq.  (1.13)  and  from  equations  similar  to 
(1.13),  [22].  Equations  (2.4)  provide  the  eigenvalues  p^,  (L  =  1,  2, 


3)  and  the  associated  eigenvectors  U.  .  and  T.,  ,.  For  Composite  W, 

1 ,  L  t  J  I  i> 

all  three  eigenvalues  p^  are  purely  imaginary  for  any  angle-ply  8, 

[22],  while  for  Composite  T  two  of  the  three  eigenvalues  are  complex  for 
|0(  less  than  71.5377°.  By  substituting  Eqs .  (2.10,  2.11)  into  the 
stress  free  boundary  conditions  and  the  interface  continuity  conditions 
Eqs.  (2.12,  2.13),  one  obtains  a  system  of  12  linear  homogeneous 
equations  for  the  constants  A.,  B. ,  A.',  ' ,  Eq.  (2.14).  The 

Ij  xj  Jj 

roots  of  the  determinant  of  this  system  are  found  by  the  method 
described  in  the  previous  section.  An  area  bounded  by  (-1  <  5  <  3,  0  < 

Ti  <  3)  was  checked.  The  roots  found  for  various  (0/0')  combinations  are 
listed  in  Tables  4  and  5.  Double  precision  was  used  in  the  calculations 
but  we  have  rounded  the  roots  in  Tables  4  and  5  to  four  digits.  Both 
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complex  and  real  roots  were  found.  Since  complex  conjugates  of  these 
values  are  also  roots  of  it  was  not  necessary  to  search  the  n  < 


0  area. 

Most  interesting  of  these  roots  are  the  positive  integer  values  of 
6,  which  seem  to  consistently  appear  for  all  (0/6')  combinations  for 
both  composites.  There  appears  to  be  a  negative  real  root  for  6,  but 
there  are  no  other  complex  roots  with  a  negative  real  part.  Since  the 
negative  6  is  the  one  contributing  to  the  singular  stress,  we  list  in 
Tables  2  and  3  the  negative  6  for  various  combinations  of  (6/6')  angles. 
Also,  we  present  in  Figs.  8  and  9  the  negative  6  for  all  possible 
combinations  of  (0/0')  angles.  Curves  of  constant  6  are  given  only  in 
one  quarter  of  the  (0,  0')  plane  since  the  curves  in  the  remaining  three 
quarters  are  a  repetition  of  the  curves  shown.  The  negative  6  values 
also  appear  to  be  simple  roots  of  Eq.  (2.20),  and  hence  a  of  Eq.  (2.19) 

•it 

is  unique  up  to  a  multiplicative  constant,  say  k  .  By  substituting  a 
of  (2.19)  into  (2.18),  we  may  write  Eq.  (2.18)  as 


o  .  . 


*  6  ^ 
k  r  o  .  . 
ij 


(2.27) 


where  0.  .  depends  on  0 .  It  should  be  pointed  out  that  if  6  <  0  is 


is  a  double  root  of  Eq.  (2.15)  one  would  have,  besides  the  r  singularity, 

a  singularity  of  the  form  r^(ln  r) ,  [11]. 

The  analysis  presented  here  provides  the  order  of  singularity  6 


and  .  k  ,  which  may  be  identified  with  the  stress  intensity  factor 


if  elements  of  are  normalized,  can  be  determined  only  by  solving 


the  global  boundary  value  problem.  For  instance ^  one  may  use  a  finite 
element  scheme  in  which  a  special  element,  whose  stress  is  given  by  Eq. 


■A* 

(2.27),  is  introduced  at  the  free-edge.  If  k  associated  with  a  free- 
edge  point  happens  to  be  zero  after  solving  the  global  boundary  value 
problem,  there  is  no  singularity  at  that  particular  free-edge  point. 
Therefore,  a  singularity  at  the  free-edge  point  in  plane  strain 
deformation  is  not  certain  until  the  global  problem  is  solved. 

TABLE  1 

Values  at  the  Nodes  of  a  Square  for  Composite  W,  (60/-60) 

Node  +  irijj  =  6  u^  +  iv^^  = 

1  2.94158674804  +  1.775502774707  i  204801.982842  +  47335.881916  i 

2.94158674902  +  1.775502774707  i  -64582.9527519  -  285394.704755  i 

2.94158674902  +  1.775502774804  i  268147.7157535  -  554779.491707  i 

2.94158674804  +  1.775502774804  i  537532.4431132  -  222048.858434  i 


TABLE  2 


Negative  Real  Roots  6  for  Terms  at  the  Free-Edge  in  Composite  W 


e 

6 

0'  =  0 

e’  =  90® 

B'  =  -0 

o** 

-3.3388  X  10*^ 

15® 

-1.3528  X  lO'^ 

-3.2814  X  10*^ 

-6.4322  X  lO'^ 

O 

o 

-2.6286  X  lO'^ 

-2.8682  X  10*^ 

-1.1658  X  lO'^ 

45® 

-9.6461  X  lO'^ 

-2.0575  X  10*^ 

-2.5575  X  lO'^ 

60® 

-1.9866  X  10‘^ 

-1.0519  X  10*^ 

-2.3346  X  10*^ 

75® 

-2.9388  X  lO'^ 

-2.6785  X  lO’^ 

-8.9444  X  lO"^ 

90® 

-3.3388  X  10‘^ 

TABLE  3 


Negative  Real  Roots  6  for  r  Terms  at  the  Free-Edge  in  Composite  T 


0°  .  -5.4148  X  10 


15®  -5.5587  X  lO"^  -5.2192  x  lo'^  -2.5363  x  10 
30®  -5.8892  X  lo'^  -4.4295  x  lo"^  -2.2505  x  10 
45®  -1.8423  X  10“^  -3.0453  x  lo'^  -3.8593  x  10 
60°  -3.4756  X  lo"^  -1.4602  x  lO'^  -3.1271  x  10 
75®  -4.8646  x  lO"^  -3.4296  x  lO"^  -1.1217  x  10 
90®  -5.4148  X  lO'^  .  . 


Complex  Roots*  6  for  Terms  at  the  Free-Edge  in  Composite 


Complex  Roots*  6  for  r  Terms  at  the  Free-Edge  In  Composite 


computer  subroutine  that  checks  a  square  for  possible  complex  roots 


Chapter  III 


THREE-DIMENSIONAL  DEFORMATION 

з. 1  INTRODUCTION 

The  analyses  presented  in  Chapter  II  assumed  that  the  displacement 

и. ,  and  hence  the  strain  e..  and  the  stress  0.,,  are  functions  of 

1  ij  ij 

and  only.  In  other  words,  we  assumed  a  two-dimensional  plane 

strain  deformation.  If  the  deformation  is  three-dimensional,  u.  would 

1 

depend  on  as  well.  Before  we  derive  a  first  order  approximation 

for  three-dimensional  deformations,  we  point  out  that  the  analyses 
presented  in  Chapter  II  remain  valid  for  certain  classes  of 
displacements  which  depend  on  x^ . 

Firstly,  we  may  superimpose  the  displacement  given  in  Eq.  (2.1)  by 

u.  =  U.  +  «. .X.  (3.1 

1  1  ij  J 

where  U.  and  u.,  are  constants  and  u..  is  an  antisymmetric  matrix. 

1  ij  ij 

The  displacement  given  by  Eq.  (3.1)  represents  a  rigid  body  translation 
and  rotation,  and  hence  contributes  nothing  to  the  strain  and  stress. 
Thus,  even  though  x^  is  present  in  Eq.  (3.1),  the  analyses  in  Chapter 

II  remain  valid  when  Eq.  (3.1)  is  superimposed  on  Eq.  (2.1). 

Secondly,  we  will  show  that  the  homogeneous  strain  produced  by  the 
displacement  field 
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where  and  c.^  are  constants  can  also  be  produced  by  the 

formulation  in  Chapter  II  by  letting  6=0.  Indeed,  when  we  put  6=0 
in  Eq.  (2.6),  u.  is  linear  in  x,  and  x.,  and  hence  t..  are 

constants.  To  insure  that  assume  real  values  with  6=0,  we  let 

Bt  =  A,.  ,  (L  =  1,  2,  3).  Since  are  complex,  we  have  six  real 
Ij  L  ij 

arbitrary  constants.  By  introducing  another  six  real  arbitrary 

constants  a^,  b^,  c^,  (i  =  1,  2),  we  may  rewrite  Eq.  (2.6)  for 

6  =  0  as 


’^1  ^2^2 


“2  =  ^  V2 


'^3  =  ‘=2*2 


(3.3) 


If  a^,  b^,  c^  are  all  non-zero,  Eq.  (3.3)  provides  all  six  strain 


components  except  £23'  Now,  consider  the  following  special  case  of 


Eq.  (3.3); 

Ui  =  U2  =  0,  u^  =  CjXj  +  0^X2  (3.4) 

It  is  not  difficult  to  show  that  the  strain  obtained  from  Eq.  (3.2)  and 
Eq.  (3.4)  is  identical.  Thus  the  deformation,  Eq.  (3.2),  which  depends 


on  x^  is  a  special  case  of  the  deformations  considered  in  Chapter  II. 
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3.2  FIRST  ORDER  EXPANSION 


Let  the  displacement  near  the  origin  in  Fig.  2  be  a  function 
of  and  x^.  For  points  near  the  x^  =  0  plane,  we  expand 

in  the  following  series 


11  1  3  1  3 


(3.5) 


where  ,  (k  =  0 ,  1 ,  2 ,  .  .  . )  are  functions  of  and 

only.  The  strain  is,  by  Eq.  (1.1) 

+  .  .  .  (3.6) 

where  6^^  is  the  Kronecker  delta.  We  assume  that  the  displacement 

(k) 

u.  is  continuous  and  bounded.  Then  u.  ^  are  continuous  and 
1  1 

bounded.  In  Eq.  (3.6),  the  only  terras  which  may  contribute  to  a 

singularity  in  must  therefore  come  from  the  derivatives  of 

(k)  (k)  fkl 

u.  ,  i.e.  from  the  terms  u.  .  and  u .  .  .  We  also 

1  i.J  J.i 

assume  that  the  order  of  singularity  6  at  x^  =  x^  =  0,  if  it  exists, 
is  continuous  in  x^ •  This  means  that  the  order  of  singularity  of  the 
(k) 

term  u.  ,  (k  =  1,  2  .  .  .)  cannot  be  stronger  than  the  order  of 

1  >  J 

the  singularity  of  Therefore,  to  the  first  order  of 

approximation  as  (x^,  x^,  x^)  -♦  (0,  0,  0)  the  strains  are 


^-*1 
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The  terms  in  the  first  parentheses  provide  singular  (infinite)  strains 
while  the  terms  in  the  second  parantheses  yield  a  finite  value  because 


u.  is  continuous  and  bounded.  Let 
1 


=  c.  +  .  .  .  (3.8) 

where  c^,  (i  =  1,  2,  3)  are  constants  and  the  dots  represent  terms 
which  vanish  when  (x^,  x^)  (0,  0).  Hence,  to  the  first  order  of 

approximation,  we  may  write  Eqs .  (3.5)  as 

^i  =  +  c.x^  +  .  .  .  (3.9) 

The  term  is  the  two-dimensional  deformation  and  is  identical 

1 

to  the  right  hand  side  of  Eq.  (2.1)  while  the  term  is  the  first 

order  approximation.  However,  as  we  discussed  in  Eq.  (3.2),  the 
deformation  for  =  c^x^,  u^  =  c^x^  can  be  included  in  the 

two-dimensional  deformations.  We  may  therefore  omit  c^  and  c^ 

without  loss  of  generality,  and  write  Eq.  (3.9)  as 

Ui  =  u^fCz)  +  ^26^3X2  (3.10) 

where  =  c^  is  the  strain  in  the  x^-  direction. 

Thus,  to  the  first  order  of  approximation,  the  stress  singularity 
at  a  wedge  corner  in  three-dimensional  deformations  can  be  analyzed  by 
considering  the  two-dimensional  plane  strain  problem  superimposed  by  a 
uniform  extension  in  the  x^-direction .  In  the  next  chapter. 
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Chapter  IV,  we  will  consider  stress  singularities  at  a  free-edge  point. 
In  Chapter  V  the  singularities  at  a  contact-edge  point  will  be 
discussed. 

The  seemingly  innocent  appearance  of  the  uniform  extension  term,  as 
we  will  see  later,  makes  the  stress  at  the  free-edge  point  inherently 
singular  for  certain  composites.  Moreover,  the  singularity  is 
logarithmic . 


Chapter  IV 


THREE-DIMENSIONAL  FREE-EDGE  SINGULARITY  ANALYSIS 
The  problem  studied  here  is  the  same  as  shown  in  Fig.  2  except  that 
now  uniform  extension  is  added;  thus  making  this  problem  a  first 

order  approximation  in  three  dimensions.  We  will  first  assume  that  a 
homogeneous  (uniform)  stress  solution  exists  for  a  specified  t^.  Two 

methods  will  be  used  to  analyze  this  problem.  They  give  the  same 
solution  for  the  cases  of  cross-ply  and  (6/-0)  composites.  For 
composites  which  are  neither  cross-ply  nor  (0/-0),  a  uniform  stress 

* 

solution  due  to  the  extension  does  not  exist.  For  these  composites, 

we  find  a  solution  in  which  (In  r)  terms  appear  in  the  expressions  for 
stress  and  displacement.  These  will,  of  course,  lead  to  singularities 

at  r  =  0.  Singularities  of  the  type  r^  (-1  <  Re(6)  <  0)  found  in 

the  analysis  of  the  two  dimensional  plane  strain  problem,  Chapter  II, 

may  still  exist,  and  should  be  superimposed  on  the  solution  obtained  here. 


4.1  HOMOGENEOUS  STRESS  SOLUTION:  METHOD  1 

In  this  section  we  assume  that  a  homogeneous  stress  solution  exists 

due  to  a  uniform  extension  in  the  x_-direction.  Hence  the  stress  o..  and 

3  ij 

the  strain  are  independent  of  x^  and  x^ •  By  using  the  notations 
discussed  in  Section  1.2  along  with  Eq.  (1.8,  1.9,  1.10)  one  can  obtain 


4r 

For  an  exception  see  page  45 
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the  relations  between  the  strain  £..  and  the  stress  o  .  in  the 

IJ  IJ 

following  form 

'l  '  *11“!  *  ‘n”!  *  "13'’3  +  =14°4  *  “15“5  *  =‘l6“6  “ 

‘3  “  “si"!  +  "32'’2  *  =33'’3  *  ^3«“4  *  =3S°5  *  '‘36'’6 

‘5  =  *51“!  *  *52'’2  *  "53“3  *  *54°4  *  *S5°5  +  "56°6  “ 

where  the  expressions  for  needed  in  the 

following  analysis,  and  hence,  are  omitted  here.  Since  the  material  is 
symmetric  with  respect  to  the  (x^,  x^)  plane,  [25,  26) 

^14  "  ^16  ^  ®24  "  ^26  "  ®34  “  ^36  "  ®54  "  ®56  ^  ° 

By  applying  the  above  material  symmetry  properties,  along  with  the 
stress  free  boundary  conditions  at  (>  =  ±ii/2,  Eq.  (2.12),  one  obtains 


"l  =  "12®2  ^  ^13°3 

^3  =  ^32^2  ^  ^33°3 

*^5  ~  ^52°2  *  ^53^3 

Solving  for  from  Eq.  (4.6),  and  eliminating  in  Eqs.  (4.5, 
4.7)  we  have 

‘’3 

^5 


(4.5) 

(4.6) 

(4.7) 


®3'^®33  '  ^®32'^^33^°2 

(4.8) 

~  ^3®13^®33  ^®12  '  ®13®32^®33^®2 

(4.9) 

■  ^®52  ’  ®53®32^®33^°2 

(4.10) 

By  setting 


^13^^33 


^12  “  ®12  ■  ®13®32^®33 


*^5  ~  ®53'^®33 


^52  *  ®52  “  ®53®32^®33 


(4.11) 
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I 


Eq.  (4.9)  and  Eq.  (4.10)  can  be  rewriten  as 


^12°2 


^5  =  ^5  =  3  ®52°2 


(4.12) 

(4.13) 


The  interface  continuity  conditions  of  Eq.  (2.13a)  at  0  =  0  are 
equivalent  to 


(Ell  =  [E3I  =  (£5)  =  0  (4.14) 

Using  Eqs.  (2.13b)  and  (4.14)  Eqs.  (4.12,  4.13)  become,  respectively, 

[Rl2la2  t^ll^  =  °  (^-15) 

=  0  (4.16) 

For  cross-ply  composites,  i.e.  (0/90)  composites,  = 

R^^'  ~  0  because  we  have  assumed  that  material  1  is  identical  to  material  2. 
Hence  Eq.  (4.16)  is  automatically  satisfied  and  Eq.  (4.15)  yields 

O2  =  -^3(^11/1^12! 

For  (0/-0)  composites,  Rj2  ~  ^12^  ^1  ~  ^l'’  Equation  (4.15) 

is  automatically  satisfied  and  Eq.  (4.16)  provides  a^: 

®2  =  -^(^51/(^521 

For  other  (0/0')  combinations,  Eqs.  (4.15)  and  (4.16)  contradict  each 
* 

other.  This  indicates  that  the  assumed  homogeneous  stress  solution 

due  to  a  prescribed  extension  does  not  exist  for  other  (6/9')  composites. 

When  a  homogeneous  solution  exists,  is  obtained  from  Eq.  (4.17) 

or  (4.18),  O3  is  from  Eq.  (4.8),  while  is  arbitrary.  If  we 

denote  by  02^^^’  *^3^^^  values  of  and  O3  obtained  from 


* 

For  an  exception  see  page  45 
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Eqs.  (4.17  or  4.18)  and  (4.8)  we  may  write  the  homogeneous  solution  as 


o .  . 

ij 


0  0  0 

1 - 

0 

0 

0 

0  0 

^  “1 

0  0  1 

0  0 

0  1  0 

(4.19) 


where  is  an  arbitrary  constant. 


A  similar  equation  applies  to 


4.2  HOMOGENEOUS  STRESS  SOLUTION:  METHOD  2 

The  preceeding  approach  fails  to  give  a  solution  when  Eqs.  (4.15) 
and  (4.16)  are  incompatible.  The  following  method,  which  is  a  modified 
version  of  the  method  in  Chapter  II,  provides  a  means  to  solve  the  cases 
for  which  a  homogeneous  stress  solution  does  not  exist.  However,  we 
will  first  apply  the  method  to  find  the  homogeneous  stress  solution  when 
it  exists,  and  see  if  this  method  produces  the  same  solution  as  in  the 
previous  section. 

To  take  into  account  the  uniform  extension  in  the  x^  direction, 


Eq.  (2.1)  is  replaced  by  Eq.  (3.10)  so  that 


u.  =  u.f(Z)  +  6.3E3X3 


Of.  =  X  .  .df(Z)/dZ  c.  .33^3 


By  Eqs.  (1.1,  1.2)  one  obtains  expressions  for  p,  x..,  u, 

1 J  K 


(4.20) 

(4.21) 


which  are 


the  same  as  those  for  the  two  dimensiunal  plane  strain  problem, 

Eq.  (2.4).  cihoosing  f(Z)  to  have  the  form  defined  in  Eq.  (2.5),  Eqs. 
(4.20,  4.21)  may  be  rewritten  as 


'  *  Vij.A*’  *  ^iJ33‘3 


o .  . 
IJ 


(4.23) 
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where  A.  and  B.  are  complex  constants.  Equations  (4.22,  4.23) 

ii  ij 

differ  from  Eqs.  (2.6,  2.7)  in  the  last  term.  Since  Eqs .  (2.6,  2.7)  can 
be  written  in  the  form  of  Eqs.  (2.17,  2.18)  when  6  is  real,  we  rewrite 
Eqs.  (4.22,  4.23)  for  real  6  as 

u.  =  r'^^  Z  +  aj^Im(u.  }/(l  +  6) 


"3^i3’'3 


(4.24) 


^  "ij33=3 


(4.25) 


where  a^  and  a,,  are  real  coefficients.  By  applying  the  stress  free 

Ij  Li 

boundary  conditions  specified  in  Eq.  (2.12)  at  ^  =  ±11/2,  and  the 
interface  continuity  condition  at  0  =  0,  Eq.  (2.13),  one  obtains  12 
linear  equations  for  a^^,  a^,  which  can  be  written  as 


r  K(6)a  =  £3^ 

where  b  is  a  column  matrix  whose  elements  consist  of  c..„„  and 

ij33 

^'ij33‘  special  case  when  -  0,  Eq.  (4.26)  reduces  to 

Eq.  (2.19). 

Since  the  right  hand  side  of  Eq.  (4.26)  is  independent  of  r, 
Eq.  (4.26)  holds  only  for  6=0: 


(4.26) 


K(0)a  =  E,b 


«%»  <»V 


(4.27) 

Assuming  that  (^.27')  has  a  solution  for  a  whose  elements  are  a^  ,  a., 

a.',  al ' ,  and  substitute  the  results  into  Eqs.  (4.24,  4.25)  with  6  = 

L  L 


0  we  obtain 


u.  =  X  {aLRe(u^^j^ZL)  +  aj^Im(u.  +  E36.3X3 

Oi.  =  Z  '=ij33^3 


(4.28) 

(4.29) 
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where  Eq.  (2.8)  has  been  used  in  obtaining  Eq-  (4.28).  Equation  (4.29) 
shows  that  is  homogeneous. 

As  we  pointed  out  in  Section  2.1,  the  determinant  of  K(6)  vanishes 
for  6=0.  Hence,  a  solution  to  Eq.  (4.27)  exists  if  and  only  if,  [27] 


T 

4  b  =  0 


(4.30) 


.T 


where  I  is  a  left  eigenvector  of  K(0) : 

4^K(0)  =  0  (4.31) 

It  turns  out  that  there  are  two  left  eigenvectors  of  K(0). 

Equation  (4.30)  then,  must  be  satisfied  for  both  1.  For  the  (0/90)  and 
(0/-0)  composites,  numerical  solutions  indicated  that  Eq.  (4.30)  holds, 

and  Eq.  (4.27)  has  a  particular  solution  ^wo  arbitrary 

solutions  a^^^  and  a^^^.  Hence 


a  =  633^^^  +  (4.32) 

where  and  are  arbitrary  constants.  Substitution  of  Eq.  (4.32) 
into  Eqs .  (4.28,  4.29)  yields 


1  3  1  1  1  2  1 


0..  = 

ij  3  ij  1  ij 


(P)  _ 


ij 


=  Z  {a/P^Re(T..  ,)  +  i/P^Im(T..  .)>  +  c..,., 
L  iJ,L^  L  ij33 


(4.33) 

(4.34) 


where 

u  (P)  = 

i 

£  *  ‘i3-3 

(4.35a) 

(n) 

u.  = 

1 

Z  taL^"^Re(u^^j^Z^)  +  aj^^"hm(u.  =  1.  2) 

(4.35b) 

(4.35c) 


(4.35d) 
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(2) 

The  reason  a..  is  absent  in  Eq.  (4.34)  is  because'we  have  chosen 
ij 

and  a^^^  such  that  is  a  rigid  body  rotation,  and 

(2)  (2) 
hence  o..  associated  with  u.  is  zero, 
ij  1 

Numerical  solutions  for  the  composites  considered  here  show  that 
Eq.  (4.34)  is  identical  to  Eq.  (4.19)  for  (0/90)  and  (0/-0)  composites. 
When  the  composite  is  neither  (0/90)  nor  (9/-0)  Eq.  (4.30)  does  not 

.J. 

hold,  and  there  is  no  solution  for  a  from  Eq.  (4.27).  In  the  next 


section  we  will  treat  this  case  by  modifying  the  approach  of  this  section. 


4.3  LOGARITHMIC  SINGULARITY 

If  no  solution  exists  for  Eq.  (4.27)  the  assumption  of  uniform 
stress  due  to  a  uniform  extension  is  appparently  invalid.  In  such 


a  case,  instead  of  using  Eqs.  (4.22,  4.23)  we  use  the  following  solution: 


O..  =  ^  {  I  [Ajt.  B.T..  Z^^]}  +C..E 

IJ  ‘  L  1J,L  L  L  ij,L  L  '  ij33  3 


(4.37) 


where  Aj^,  are  now  functions  of  6.  It  can  be  shown  that  Eqs. 

(4.36,  4.37)  satisfy  Eqs.  (1.1,  1.2,  1.4)  with  u.  ^  and  t.,  given 

1 ,  L  ij  1 L 

by  Eq.  (2.4).  Before  we  substitute  Eqs.  (4.36,  4.37)  into  the  boundary 
and  interface  conditions,  we  apply  the  equivalence  between  Eqs.  (4.22, 
4.23)  and  Eqs.  (4.24,  4.25)  to  Eqs.  (4.36,  4.37).  The  result  is 


*  ‘3*13*3 


(4.38) 


For  an  exception  see  page  45 
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where  it  is  understood  that  all  quantities  on  the  left  hand  side  of 
Eqs.  (4.45,  4.46)  are  evaluated  at  6  =  0.  Equations  (4.45,  4.46) 
consist  of  24  equations  for  g  and  da/d6.  If  a  solution  exists, 
substitution  of  g  and  dg/d6  back  into  Eq.  (4.41)  with  6=0  provides  the 


desired  solution. 
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Before  we  discuss  the  solution  of  Eqs.  (4.45,  4.46)  in  the  next 
section,  we  write  Eqs.  (4.40)  and  (4.41)  in  full  with  6=0: 
u.  =  (In  r)Z  1 

+  (daj^/d6)Re(u^  ^Z^^)  +  (drj_/d6)Im(u .  (^-^7) 

0..  =  (In  r)  I  [aLRe(t..  l)  + 

+  I  [a^Re(T..^^ln4^)  +  tj^Im(t .  .  ^  ^InC^^) 

+  (da^/d6)Re(T  +  (da^/d6)In)(T  .  ,  ^^)]  +  (4.48) 

We  see  that  has  a  logarithmic  singularity.  Again,  Eq.  (2.8)  has 
been  used  in  obtaining  Eq.  (4.47). 

4.4  SOLUTION  FOR  STRESS  AND  DISPLACEMENT 

The  system  of  Equations  (4.45,  4.46)  has  a  unique  solution  for  a 

'V 

if,  (see  [10]) 

d^liK||/d6^  ^  0,  N  =  n  -  m,  (4.49) 

where  n  and  m  are,  respectively,  the  order  and  rank  of  K.  For  the 
composites  considered  here,  N  =  2.  However,  it  is  rather  difficult  to 
prove  or  disprove  Eq.  (4.49)  analytically  or  numerically  in  view  of  the 
fact  that  K  is  a  12  X  12  matrix.  Instead,  we  regard  Eqs.  (4.45,  4.46) 
as  a  system  of  24  equations  for  ^  and  d_a/d6,  and  solve  the  system 
numerically.  We  find  numerically  that  a  is  unique  while  da/d6  has  a 
particular  solution  and  two  arbitrary  solutions. 

Noting  that  N  =  2,  one  can  see  that  K  has  two  right  eigenvectors 

a^^^  and  a^^^  such  that 

•'v*  /V 

Ka^"^  =0,  (n  =  1,  2)  (4.50) 
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If  a  is  the  unique  solution  of  both  Eqs .  (4.45)  and  (4.46),  it  mu^t  also 
be  a  solution  of  Eq.  (4.45),  and  hence  a  is  proportional  to  a  right 

eigenvector.  Without  loss  of  generality,  let  be  the  eigenvector 

to  which  a  is  proportional,  i.e., 

a=kS3a^l>  (4.51) 

where  k  is  uniquely  determined  if  a^^^  is  properly  normalized.  The 
fact  that  da/d6  has  two  arbitrary  solutions  is  obvious  from  Eq.  (4.46) 
because  the  coefficient  of  da/d6  is  K  which  is  singular  of  order  2.  If 

£3(da/d6)^^^  is  a  particular  solution  of  da/d6,  we  have 


da/d6  =  E,(da/d5)^^^  +  o  a^^^  +  a  a^^^  (4.52) 

where  Oj^,  are  arbitrary  constants.  With  Eqs.  (4.51,  4.52),  Eqs. 
(4.47,  4.48)  can  be  rewritten  as 

u.  =  k£,[(ln  r)u.^^^  +  u.^^^1  +  +  o.u.^^^  +  (4.53) 

X  X  xx  ^x 

o..  =  k.3(<l„  *  «,.<«))  *  ,3  <,..<■»  +  (4,54) 

where  u.^^\  u.^^\  and  are  identical  to  the  ones  defined 

11  ij 

in  Eqs.  (4.35b,  4.35d),  while 

u/P^  =  E  {(da^/d6)^P^Re(u.  ^Z.)  +  (  dal/dS)  ^P^Im(u .  .Z.)} 

X  Ij  Lt 


+  «13X3  (4.554) 

u.<*>  =  I  -  1)] 

-V  (  11 

+  a^^^  •'lm(u.  ^^Zj^dn  *  DD  (4.S5b) 

J  {(daj^/d6)^P^Re(x.^  j^)  +  (daj^/d6)  ^P^Im(T  j^)} 

*  ‘^ij33 


C4.55C) 
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”ij“’  ■  ^  V  *  «L»  <'>■”'1) 


(2) 

Again  the  reason  is  missing  in  Eq.  (4.54)  is  due  to  the  fact 


that  is  a  rigid  body  rotation. 

Although  the  solution  obtained  here  is  for  composites  which  are 
neither  (0/90)  nor  (0/-0),  application  of  the  present  solution  to  (0/90) 
or  (0/-0)  composites  yields  k  =  0.  Hence  a.  =  £  by  Eq.  (4.51),  and  the 
solution  for  da/d6  from  Eq.  (4.46)  is  identical  to  the  solution  for  a  in 

Eq.  (4.27).  It  follows  that  (da/d6)^^^  of  Eq.  (4.52)  is  identical  to 

a^^^  of  Eq.  (4.32),  and  that  u.^^\  o . . in  Eqs.  (4.55a,  4.55c) 

-v  1  ij 

and  Eqs.  (4.35a,  4.35c)  are  also  identical.  Thus,  the  solution  obtained 
in  Eqs.  (4.53,  4.54)  reduces  to  that  given  in  Eqs.  (4.33,  4.34)  when  the 
composite  is  (0/90)  or  (0/-0). 


For  composites  which  are  neither  (0/90)  nor  (0/“0),  k  0.^  We 

see  from  Eq.  (4.54)  that  o.,  has  a  logarithmic  singularity  unless  t 

1 J  j 

=  0.  Therefore,  unless  happens  to  be  zero  in  a  three-dimensional 


deformation,  the  stress  is  inherently  singular  for  composites  which  are 

neither  (0/90)  nor  (0/-0).^  Moreover,  the  singularity  is  logarithmic. 
Since  the  larger  the  value  of  k  the  stronger  the  logarithmic 
singularity,  k  may  be  regarded  as  the  "logarithmic  stress  intensity 

factor".  It  should  be  pointed  out  that  the  singularity  of  k  r  ,  (6 
<  0) ,  as  analyzed  in  Chapter  II,  may  still  exist  for  all  composites. 

However,  the  determination  of  the  intensity  factor  k  requires  a 
global  solution  while  k  in  Eq.  (4.54)  does  not. 


For  an  exception  see  page  45 
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If  we  replace  the  arbitrary  constant  by  another  arbitrary. 


constant  defined  by 


=  -kt^ln 


(4.56) 


Eq.  (4.54)  can  be  rewritten  as 


a  =  E3{k  ln(r/ep  +  o.  (0)} 


(4.57) 


where 


0  .  .  (<>)  =  ko  .  .  +  o  .  . 

ij  iJ  ij 


(0)^^  (P) 


(4.58) 


Thus  Pj  must  have  the  same  physical  dimension  as  r.  We  see  from  Eqs . 

(4.35c,  4.35d)  that  and  are  constants.  For  both 

ij 


composites  considered  here,  ‘  kas  the  form 


(1)  _ 


0  0  0 
0  0  1 
0  10 


(4.59) 


On  the  other  hand,  we  see  from  Eqs.  (2.9,  4.55d)  that  depends  on 

0,  (Fig.  3).  Hence  depend  on  (p.  In  Tables  6-9 

we  list  the  logarithmic  stress  intensity  factor  k  and  ^nd  Oj^j'(^) 

on  the  interface  (0  =  0),  and  on  the  free-edge  surface  (0  =  ±90°). 

Notice  that  3^^  =  =  0  at  0  =  90°.  Similarly, 

=  ^12'  ~  ^is'  “  0  at  0  =  -90°.  Hence  these  components  are  not  listed 

in  the  tables.  Also,  since  O21'  ~  °21’  °23'  ~  °23 

0  =  0,  only  a^^CO),  021'-®^’  ®23^^^  listed.  All  in  Tables  6-10 


have  the  unit  of  10  psi. 
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For  (0/90)  or  (6/-0)  composites,  k  =  0,  and  Eq.  (4,54)  may  be 


rewritten  as 


where 


0 . .  =0,0..  +  £-0 . . 

ij  1  ij  3  ij 


c  .  .  =  0  .  . 
ij  ij 


(4.60) 


(4.61) 


is  now  independent  of  Referring  to  Eq.  (4.19),  we  see  that  the  only 
non-zero  components  for  o^^  are  and  The  numerical 

•V 

calculations  of  for  (0/90)  composites  are  included  in  Tables  6-9 

while  that  for  (0/-0)  composites  are  given  in  Table  10 

We  see  from  Table  6  that  for  the  (0/15)  composite  in 

Table  6  is  many  times  larger  than  and  .  In  other  words,  the 

tensile  stress  at  the  interface  is  many  times  larger  than  the  applied 

axial  extensional  stress.  Notice  however  that  a.,  is  not  the 

ij 

(1) 

total  stress.  The  total  stress  consists  of  o . .  and  o..  as  well 
6 

as  the  k  r  and  k(ln  r)  terms. 

In  Figs.  10  and  11  we  present  the  logaritlimic  stress  intensity 
factor  k  for  all  possible  combinations  of  (0/0').  The  contour  lines  for 
constant  k  values  are  given  only  in  one  quarter  of  the  (0,  0')  plane 
since  the  contour  lines  in  the  remaining  plane  are  a  repetition  of  the 
ones  shown  in  the  figures.  As  we  pointed  out  earlier  in  the  analysis. 
Figs.  10  and  11  indicate  that  k  =  0  occurs  for  (0/90)  and  (0/-0) 
composites.  Of  course,  k  is  also  zero  for  (0/0)  composites  in  which  the 


two  materials  are  the  same  across  the  interface.  However,  we  find  an 
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additional  contour  line  for  k  =  0  in  Figs.  10  and  11  which  lies  slightly 
above  the  line  9'  =0.  This  is  unexpected  and  was  not  obvious  from  the 
analysis  presented  earlier. 


In  conclusion,  it  should  be  pointed  out  that  the  singular  stress  of 


*  6 

the  type  k  r  obtained  in  (2.27)  should  be  added  to  (4.54)  to  obtain 
the  singularities  at  the  free  edge.  Using  (4.58),  we  may  write  the 
stress  singularity  near  the  free  edge  in  the  form 

0..  =  k  r^o .  .  +  [ke(ln  r)  +  +  t-O..  (4.62) 

ij  ij  ‘  3  1'  ij  3  ij 


*  ( 1 )  * 

where  6  <  0,  k  ,  k,  a,  and  a.,  are  constants  while  a.,  and 

1  ij  ij 

depend  on  The  analysis  presented  here  provides  every  term 


on  the  right  of  (4.62)  except  k  and  which  have  to  be  determined 


by  solving  the  complete  boundary  value  problem. 

Finally,  it  should  also  be  pointed  out  that  even  though  the  layers 
are  assumed  to  be  of  the  same  orthotropic  material  for  the  numerical 
calculations  in  this  report,  the  theory  presented  here  applies  to  any 
anisotropic  layered  composite. 
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TABLE  6 

of  Eq.  (4.58)  for  Composite  W,  (0/0') 


0' 

15° 

30° 

45° 

60° 

75° 

90° 

k 

0.0022 

0.0033 

0.0042 

0.0040 

0.0025 

0 

a'li(O) 

0.0221 

0.0145 

0.0102 

0.0056 

0.0016 

0 

o'll'(O) 

-0.0222 

-0.0149 

-0.0112 

-0.0067 

-0.0021 

0 

.  . 

^22(^0) 

190.43 

29.599 

8.0106 

1.2162 

-1.3911 

-2.1000  , 

022^0) 

190.45 

29.613 

8.0208 

1.2218 

-1.3895 

-2.1000 

190.47 

29.628 

8.0310 

1.2274 

-1.3880 

-2.1000 

033(90) 

59.990 

26.216 

21.682 

20.256 

19.708 

19.559 

0-33(0) 

60.000 

26.222 

21.687 

20.258 

19.708 

19.559 

O33'(0) 

35.132 

7.1872 

3.4312 

2.2440 

1.7848 

1.6590  i 

1 

1 

033'(*90) 

35.137 

7.1903 

3.4340 

2.2461 

1.7856 

1.6590  1 

0^3(90) 

0 

0 

0 

0 

0 

: 

0  ’ 

i 

923(0) 

0 

0 

0 

0 

0 

: 

0  j 

023'(-90) 

0.0003 

0.0006 

0.0007 

0.0004 

0 

0  1 
_ i 

5,3(0) 

-0.0035 

-0.0052 

-0.0065 

-0.0063 

-0.0040 

0 

5i3'(0) 

0.0034 

0.0048 

0.0059 

0.0059 

0.0039 

0 

Oi2(0) 

0.0140 

0.0092 

0.0065 

0.0035 

_ 

0.0010 

0 
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TABLE  7 

of  Eq.  (4.58)  for  Composite  W,  (90/6') 


0' 

0“ 

15“ 

30“ 

45“ 

60“ 

75“ 

k 

0 

0.6874 

0.5030 

0.2954 

0.1297 

0.0239 

o^l(O) 

0 

-0.3032 

-1.0048 

-1.3511 

-1.1657 

-0.4910 

3ii’(0) 

0 

0.2058 

0.7794 

1.1633 

1.0875 

0.4823 

522(90) 

-2.1000 

-4.5417 

-4.9455 

-4.2661 

-3.3602 

-2.4156 

522(0) 

-2.1000 

-4.6419 

-5.2778 

-4.7129 

-3.7456 

-2.5788 

522* (-00) 

-2.1000 

-4.7473 

-5.6151 

. 

-5.1610 

-4.1303 

-2.7409 

833(90) 

1.6590 

1.1463 

1.0615 

1.2041 

1.3944 

1.5925 

033(0) 

1.6590 

1.1185 

0.9695 

1.0805 

1.2877 

1.5476 

333' (0) 

19.559 

11.311 

5.5014 

2.9133 

1.8095 

1.5966 

033' (-90) 

19.559 

7.8189 

2.9701 

1.7004 

1.3932 

1.5244 

O23C9O) 

0 

0 

0 

0 

0 

0 

523(0) 

0 

0 

0 

0 

0 

0  1 

1 

023’(-00) 

0 

0.2870 

0.3950 

0.2754 

0.1242 

0.0222 

5i3(0) 

0 

-1.0797 

-0.7902 

-0.4641 

-0.2038 

-0.3762 

1 

3i3'(0) 

0 

1.6460 

1.8195 

1.3607 

0.6910 

i 

0.1382 

3i2(0) 

0 

-0.0942 

-0.3122 

-0.4197 

-0.3621 

-0.1525 

TABLE  8 


of  Eq.  (4.58)  for  Composite  T,  (0/0') 


0' 

15“ 

30“ 

45“ 

60“ 

75“ 

90“ 

0.0349 

0.0610 

0.0729 

0.0658 

0.0394 

0 

3„(0) 

0.2203 

0.1919 

0.1335 

0.0695 

0.0192 

0 

o'li'CO) 

-0.2218 

-0.1990 

-0.1467 

-0.0817 

-0.0238 

0 

©22(90) 

62.232 

13.729 

3.3929 

0.0314 

-1.2148 

-1.5400 

o~22(0) 

67.452 

13.921 

3.5264 

0.1009 

-1.1956 

-1.5400 

©22' (-^0) 

67.672 

14.115 

3.6623 

0.1729 

-1.1753 

-1.5400 

333(90) 

40.825 

25.844 

22.950 

22.009 

21.660 

21.569 

333(0) 

40.948 

25.952 

23.025 

22.048 

21.671 

21.569 

333' (0) 

22.241 

6.0686 

2.8463 

1.7151 

1.2441 

1.1088 

©33' (-90) 

22.260 

6.0653 

2.8494 

1.7218 

1.2474 

1.1088 

323(90) 

0 

0 

0 

0 

0 

0 

523(0) 

0 

0 

0 

0 

0 

0 

©23' (-90) 

0.0036 

0.0115 

0.0118 

0.0062 

0.0011 

0 

3i3(0) 

-0.0549 

-0.0959 

-0.1145 

-0.1034 

-0.0619 

0 

©,3’(0) 

0.0498 

0.0893 

0.1076 

0.0997 

0.0612 

0 

'©12(0) 

0 . 1444 

0.1258 

0.0875 

0.0455 

0.0126 

0 

50 


0' 

0- 

j  15“ 

30“ 

022(0) 


022'(-50) 


a„(90 


0.7217 


-0.4729 


0.5122 


-1.2109 


1.1088 

1.1088 

21.569 

21.569 


0.3070  j  0.9159 


■2.3081  •  -2.1663 

\ 

\ 

-2.4380  i  -2.4991 
-2.5715  j  -2.8316 


0.8937  0.9334 

0.8481  0.8165 

t 

I 

12.089  [  5.7896 

! 

7.9469  !  2.8828 


0.2993  !  0.1331  0.0252 


-1.5256  I  -1.2991  -0.5560 

I 

1.3004  1.2090  0.5460 

’  i 


-1.5927  -0.8749  -0.0430 
-2.0120  -1.2320  -0.1958 
-2.4273  -1.5858  -0.3481 


1.0940  1.2950  1.5280 
0.9467  1.1696  1.4743 
3.0108  1.8057  1.5557 


1.6049 


1.3092  1  1.4674 


1 - - 

023(90)  : 

1 

0 

■  ■  — 1 

0 

0 

0 

0 

0 

523(0) 

0 

0 

0 

0 

0 

0 

523'(-90) 

0 

0.3258 

0.4190 

0.2847 

0.1286 

0.0237 

5i3(0) 

0 

-1.1337 

-0.8046 

-0.4701 

-0.2090 

-0.0397 

3i3'(0) 

0 

1.7992 

1.9660 

1.4739 

0.7590 

0.1559 

Oi2(0) 


-0.4353  !  -0.3707 


-0.1586 
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TABLE  10 

©22  ^22  (0/"0)  Composites 


6 

15“ 

30° 

45° 

— 

60° 

75° 

Composite 

®22  ^  ®22 

191.14 

29.887 

8.2424 

1.4187 

-1.2079 

W 

°33  ^  ®33 

35.231 

7.2214 

3.4618 

2.2766 

1.8203 

Composite 

°22  ^  °22 

71.464 

16.660 

5.9075 

2.2697 

0.8352 

T 

°33  °33 

23.032 

6.5434 

3.3082 

2.2137 

1.7821 

Chapter  V 


A  CONTACT  PROBLEM  IN  COMPOSITES 


5.1  TWO-DIMENSIONAL  PROBLEM 

In  Fig.  2,  the  edge  =  0  of  the  composite  was  assumed  stress 


free.  If  a  rigid  body  is  pressed  against  the  face  =  0,  we  would 


have  the  deformation  illustrated  in  Fig.  12.  Assuming  that  the  rigid 
body  is  smooth  and  friction  free,  the  boundary  conditions  on  x^  =  0 


and  near  the  origin  are  now  given  by 


“12  '  “■ 


”13  “  “ 


Uj  =  U  +  wx^  + 


(5.1) 


where  U  and  w  are  constants,  and  the  dots  stand  for  terms  of  order  equal 

2 

to  or  higher  than  X2  •  The  first  two  terms  for  u^  can  be  included 

in  a  rigid  body  displacement  as  we  discussed  in  Chapter  III.  Therefore, 
for  the  purpose  of  finding  the  singularity  at  the  origin,  we  may  replace 
Eq.  (5.1)  by 

Uj=O3=o^=0  (5.2) 

The  eigenvalues  Pj^  and  the  associated  eigenvectors  u^,  as 

well  as  the  expressions  for  stress  and  displacement  are  the  same  as  in 
Chapter  II,  Eqs.  (2.4-2. 7).  By  applying  the  interface  continuity 
condition  Eq.  (2.13)  at  0  =  0,  and  the  new  boundary  conditions  Eq.  (5.2) 
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at  ^  =  ±ir/2  to  Eqs.  (2.10,  2.11),  we  obtain  a  new  system  of  12  linear 
homogeneous  equations  for  A  ,  B,. ,  A-  ' ,  '  which  can  be  written  as 

where  the  hat  is  used  to  distinguish  this  new  matrix  from  the  previous 
matrix,  and  q  is  again  a  column  matrix  whose  elements  are  the 

~  Vo  ~ 

complex  constants  Aj^,  Bj^,  Aj^',  B^^' .  If  6  is  real  then  Eqs. 

(2.17,  2.18)  may  be  used.  The  system  of  equations  obtained  by  applying 
the  boundary  and  interface  continuity  conditions  to  these  equations  is 

Ka  =  0  (5.4) 

where  K  is  a  real  matrix,  and  a  is  a  column  matrix  whose  elements  are 
the  real  constants  a^,  a^,  both  systems,  Eqs. 

(5.3,  5.4),  a  nontrivial  solution  exists  if  ||K^||  =  0.  Equation  (5.3) 

is  solved  by  the  method  described  in  section  2.2  for  complex  6,  while  a 
simpler  and  less  time  consuming  method  is  used  to  find  real  6  from  Eq. 
(5.4).  Values  found  for  6  by  both  methods  are  listed  in  Tables  11-14. 
The  negative  real  values  of  6,  (-1  <  6  <  0)  lead  to  a  stress  singularity 
at  r  =  0 . 

5.2  THREE-DIMENSIONAL  PROBLEM 

As  we  discussed  in  Chapter  III,  a  three-dimensional  deformation  is, 
to  the  first  order  of  approximation,  a  uniform  extension 

Xg-direction  superimposed  on  the  two-dimensional  plane  strain 

deformation.  Similar  to  the  free-edge  problem  considered  in  Chapter  IV, 
we  will  first  assume  that  a  homogeneous  stress  solution  exists  due  to  a 
uniform  extensional  strain  e^.  The  problem  will  be  analyzed  by  two 


56 


methods.  The  first  method  enables  us  to  see  easily  if  a  homogeneous 
stress  solution  exists.  The  second  method  enables  us  to  extend  the 
analysis,  if  necessary,  to  the  case  when  a  homogeneous  stress  solution 
is  impossible. 


5.3  HOMOGENEOUS  STRESS  SOLUTION:  METHOD  1 


The  boundary  condition  u^^  =  0  on  x^  =  0  of  Eq.  (5.2)  implies 


that 


au^/ax^  =  auj’/ax^  =  o 


(5.5) 


along  Xj^  =  0.  If  the  stress,  and  hence  the  strain,  is  homogeneous, 

Eq.  (5.5)  also  applies  along  x^  =  0.  Now 

=  (au2/3x^)/2,  =  (au^'/axp/Z  (5.6) 

and  since  u^  is  continuous  along  X2  =  0, 

[Egl  “  0  (5.7) 

For  a  homogeneous  stress  solution  therefore,  Eq.  (5.7)  replaces  the 
condition  u^  =  0  along  x^  =  0.  Application  of  the  material  symmetry 

properties  of  Eq.  (4.4)  and  the  boundary  conditions  of  Eq.  (5.2)  to  Eq. 
(1.3)  for  Ej,  E^,  Eg,  and  e^  results  in 


"1  =  ^3°3 


"3  =  ®3l'’l  ^  "32^2  *33^3 


'5  =  "si'’!  "52°2  "53‘’3 


E,  =  S,,0, 
6  64  4 


(5.8) 

(5.9) 

(5.10) 

(5.11) 


Solving  for  from  Eq.  (5.9),  and  eliminating  from  the 
remaining  equations,  Eqs.  (5.8,  5.10)  we  have 
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<^3  " 

^3^®33 

-  ^"3l/"33)°l 

-(532/523)02 

(5.12) 

^1  = 

^"11  • 

^31^13^^33^°! 

^®12  ‘  ^32®13^®33^°2 

^®13'^^33^'^3 

(5.13) 

"5  = 

^^51  • 

®31®53/^33^‘’l 

^  ^"52  -  ^32^53/"33)‘^2  ^ 

("53/"33)^3 

(5.14) 

Let 


^  =  "l3/"33 


^5  =  "53/^33 


^1 

=  "11 

"3l"l3'^"33 

^2 

=  "12 

"l3"32^"33 

«51 

=  "51 

"53"31^"33 

S2 

"  "52 

"53"32^"33 

(5.15) 


If  we  apply  the  interface  conditions  Eqs.(5.7)  and  (4.14)  to  Eqs.  (5.11, 
5.13,  5.14),  we  obtain 

(5.16) 


^"64'^  =  0 


^11°1  ■  ^11  °1  ’^^1^^3 


^51°1  ■  ^51  ®1  ^^52^®2  '^^5^^3 


(5.17) 

(5.18) 


Equation  (5.16)  indicates  that  if  [s^^]  *  0,  which  is  the  case 
for  (0/90)  composites,  is  arbitrary.  Otherwise  0^  =  0 .  For  the 
composites  considered  here,  s^^  =  =  0  for  any  angle  ply.  Hence 

is  arbitrary  here,  and  will  be  denoted  by  the  arbitrary  constant 

“l- 

Another  arbitrary  solution  is  obtained  from  solving  Eqs.  (5.17, 
5.18).  If  we  let  =  o^,  we  may  solve  for  and  o^'>  while 

and  o^'  are  determined  from  Eq.  (5.12).  The  result  can  be 

written  in  the  following  form: 
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«  -  r  «  r>^2) 


(5.19) 


where 


.(p)  = 


and 


0 

0 

1 

0 

0 

I 

0 

! 

0 

0 

0 

= 

0 

0 

1 

1 - 

0 

0 

0 

_ 1 

1 

0 

1 

.(2)  . 


01^2)  0 


0  1 

0  0 


(2) 


-  -CtRilRji’  -  CR5ll'u'»/«ll‘<Sl'  ■  "sAl')  '' 

Oj«)  =  -ClR^jlR^j'  .  |R52lR„’)/<Riil!5i'  -  K5i«„') 


(5.20) 


>  (5.21) 


a3^P^  =  -(1  -  =3i'’i^^^^/=33 


°3^^^  =  '^*31®!^^^  ■"  ®32^/^33 


Similar  expressions  can  be  obtained  for  o*. 

For  (0/-0)  composites,  "^1'  ®ii  ~  ^11'  ^12  ~  *^12’  ^5  ~  ~^5 '  ’ 

and  R^^  =  ”^52' '  °1^^^  ^^’21)  reduced  to 

”1“”  '  -‘‘5/''5l-  -l”'  ■=  -''52'''S1  «•“> 

0 

Moreover,  0^^^  =  0'^^^,  =  0'^^^  and  hence  0=0'.  In  Tables  15  and  16 

we  list  both  composites  of  (0/-0)  ply. 
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For  (0/8')  composites  in  Tables  17  and  18,  while 


quite  different.  Tables  19  and  20  show 


(2) 

the  results  for  (90/6')  composites;  in  this  case  '  is  nearly 
(2)  (2)  (2) 

equal  to  O^'  ,  while  '  and  are  quite  different. 

Both  3ud  ^2^^^  Tables  15-20  have  the  unit  of  10^  psi. 


(2)  „  (P) 


For  (0/90)  composites,  =  0  and  o^''  , 


of  Eq.  (5.21)  do  not  exist.  In  fact  ”  ^52*  ~  ^5  ~  ^5 '  ~  ^ 

also,  and  Eq.  (5.18)  is  trivially  satisfied.  Equation  (5.16)  and  (5.17) 
then  yield  three  arbitrary  solutions.  Letting  0^  =  0^^,  02"  “2 

and  o^'  *  «3>  the  solution  can  be  written  as 


0  =  £30^^^  +  a  +  a  +  a 


(5.23) 


where  and  are  identical  to  the  ones  defined 

•s*  ^  -w 

Eqs.  (5.20,  5.21)  except  ^nd  which  are  replaced  by 


in 


o/'”  =  -IRil/Rn.  »,«>  =  -IBijI/R,, 


(5.24) 


and  0*-^^  is 


0^(3)  0 


0^3)  =  '  0 


0 

0  0 


(3)  ! 


(5.25) 


where 


0,«>  =  Bjj'/Rii, 


„  (3)  _  „(3)^  , 

3  ~  1  31^  33 


(5.26) 
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For  o',  the  components  of  a ' and  are  given  by 

O^.CP)=0.  o^-(2>=0. 

O,  '^  =  I/S33  .  0,  =  -s,,  /s,,  , 


32  '  33 


°2  ^  ^  ‘^31  ^®33 


(5.27) 


In  Table  21  we  list  the  result  for  (0/90)  composites.  For  the 

composites  considered  here  0^  ''  =  o' 

It  appears  that  Eqs .  (5.17,  5.18)  always  have  a  solution  unless  one 
of  the  following  relations  is  satisfied: 

0  =  Rii  =  =  [R12I  *  tRil  (5.28) 

0  =  R51  =  Rg^'  =  (R52I  *  IR5I  (5.29) 

R11/R5I  =  Rji’/Rgi'  =  **  [Ril/fRsl  (5-30) 

If  anyone  or  more  of  the  above  three  equations  is  satisfied,  Eq.  (5.17) 
and/or  (5.18)  are  inconsistent  and  a  solultion  does  not  exist.  For  the 
composites  considered  here,  none  of  Eqs.  (5.28,  5.29,  5.30)  are 
satisfied  for  all  possible  combinations  of  (9/9')  angles.  Therefore,  a 
homogeneous  stress  solution  always  exists  for  a  uniform  extension  E^. 

5.4  HOMOGENEOUS  STRESS  SOLUTION:  METHOD  2 

Equations  for  stress  and  displacement  are  the  same  as  Eqs.  (4.28, 
4.29).  Upon  applying  the  boundary  conditions  of  Eq.  (5.2)  at  0  =  iir/S, 
and  the  interface  continuity  conditions  of  Eq.  (2.13)  at  ^  =  0,  we 
obtain  the  system  of  equations 

Ka  =  bE., 


(5.31) 
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which  have  a  non-trivial  solution  if  2  b  =  0,  where  ^  K  = 

T 

the  composites  considered  here  i  b  =  0  for  all  (0/0')  ply, 
numerical  solurions  agree  with  those  obtained  by  method  1. 


I 


For 

and  the 
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TABLE  11 


Negative  Real  Roots  6  for  r  Terms  at  the  Contact-Edge  in  Composite  W 


6 

6 

0'  =  0 

6’  =  90“ 

0'  =  -0 

15° 

-9.3671  X  lO'^ 

-4.6536  X  10  ^ 

-4.6679  X  10‘^ 

30“ 

-1.2238  X  lO'^ 

-6.0849  X  lO'^ 

-6.8300  X  lO'^ 

45“ 

-3.6434  X  10‘^ 

1 

-1.8759  X  lO'^ 

-1.7190  X  10‘^ 

60“ 

-5.5918  X  lO"^ 

-3.0693  X  lO'^ 

-2.0388  x  lO'^ 

75“ 

-3.8771  X  10‘^ 

-2.3686  X  lO'^ 

-1.1389  X  lO"^ 

TABLE  12 


Negative  Real  Roots  6  for  r^  Terms  at  the  Contact-Edge  in  Composite  T 


— 

6 

0 

0'  =  0 

0'  =  90“ 

0'  =  -0 

15“ 

-2.5704  X  lO'^ 

-1.1110  X  lO'^ 

-1.2894  X  lO’^ 

30“ 

-1.9977  X  lO’^ 

-8.6484  X  lO’^ 

-1.0695  X  lO'^ 

45“ 

-5.0588  X  lO'^ 

-2.3031  X  lO’^ 

-2.2038  X  lO'^ 

60“ 

-7.1066  X  lO'^ 

-3.5292  X  lO’^ 

-2.3995  X  10*^ 

75“ 

-4.6735  X  lO'^ 

-2.6708  X  lO'^ 

-1.3011  X  lO'^ 
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TABLE  13 


5 

Complex  Roots  6  for  r  Terms  at  the  Contact-Edge  in  Composite  W 


0 

«1 

— 

Or 

15/-15 

-.005 

.995  +  .100  i 

1.996  +  .120  i 

2.994  +  .180  i 

0/15 

-.001 

.999  +  .055  i 

2.0  +  .053  i 

2.998  +  .099  i 

90/15 

-.001 

.998  +  .065  i 

1.998  +  .049  i 

;  60/-60 

-.204 

.701  +  .721  i 

1.677  +  .882  i 

2.449  +  1.453  i 

0/60 

-.056 

.896  +  .413  i 

2.007  +  .359  i 

2.567  +  .919  i 

90/60 

-.037 

.840  +  .524  i 

1.985  +  .271  i 

2.657  +  1.133  i 

TABLE  14 


Complex  Roots  6 


for  r" 


Terms  at  the  Contact -Edge  in  Composite  T 


e 

«2 

o» 

15/-15 

-.013 

.989  +  .156  i 

1.994  +  .185  i 

2.993  +  .268  i 

0/15 

-.003 

.997  +  .081  i 

1.999  +  .081  i 

2.997  +  .150  i 

90/15 

-.001 

.996  +  .103  i 

1.997  +  .078  i 

60/-60 

-.240 

.649  +  .771  i 

1.621  +  .948  i 

2.359  +  1.558  i 

0/60 

-.071 

.887  +  .436  i 

2.030  +  .363  i 

90/60 

-.035 

.803  +  .575  i 

1.987  +  .274  i 

2.572  +  1.246  i 

if 

zero  and  nositive  inteeprs  are  also  roots  for  fi 


6 


TABLE  15 


^(P)^  (5.21,  5.22)  for  Composite  W,  (0/-0) 


0 

15* 

o 

o 

45° 

o 

o 

75.° 

-28.239 

-5.8883 

-1.7000 

-0.2970 

0.2543 

03<P>=,3'«» 

-3.4539 

0.8866 

1.7000 

1.9725 

2.0795 

o/2)  =a^.(2) 

0.1477 

0.1970 

0.2063 

0.2093 

0.2106 

03«)=,3’<^> 

0.2024 

0.2120 

0.2137 

0.2143 

0.2146 

TABLE  16 


of  Eqs,  (5.21,  5.22)  for  Composite  T,  (0/-0) 


0 

15° 

30° 

45° 

o> 

o 

o 

75° 

a/P)  =  dj'^P^ 

-15.948 

-4.4035 

-1.6200 

-0.6304 

-0.2332 

0  (P^  =  0  '^P^ 

3  3 

2.4089 

1.7733 

1.6200 

1.5655 

1.5436 

0^2)  =  ,^.(2) 

0.2232 

0.2643 

0.2742 

0.2778 

0.2792 

0  (2)  =0  -(2) 

3  3 

0.2886 

0.2863 

— 

0.2858 

0.2792 

0.2855 
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TABLE  17 

Composite  W,  (0/0') 


0' 

15“ 

o 

o 

45“ 

60“ 

75“ 

o/o) 

-28.206  1 

i 

-5.8551 

-1.6667 

-0.2637 

0.2877 

„^.(p) 

-28.239  ! 

1 

-5.8883 

-1.7000 

-0.2970 

0.2543 

03'P> 

14.077 

18.770 

19.650 

19.945 

20.060 

„^.(p) 

-3.4539 

0.8866 

1.7000 

1.9725 

2.0795 

o,<« 

0. 1478 

0.1971 

0.2063 

0.2094 

0.2106 

0.1477 

0.1970 

0.2063 

0.2093 

0.2106 

o,«> 

0.2410 

0.2514 

0.2533 

0.2540 

0.2542 

»3'“> 

0.2024 

0.2120 

0.2137 

0.2143 

0.2146 

L«l« 


TABLE  18 


of  Eq.  (5.21)  for  Composite  T,  (0/0') 


e' 

15“ 

30“ 

45“ 

60“ 

75“ 

o/P) 

-15.565 

-3.9220 

-1.1149 

-0.1169 

0.2837 

-15.948 

-4.4035 

-1.6200 

-0.6304 

-0.2332 

C3<P) 

17.642 

20.902 

21.688 

21.967 

22.079 

2.4089 

1.7733 

1.6200 

1.5655 

1.5436 

0.2245 

0.2660 

0.2760 

0.2796 

0.2810 

0.2232 

0.2643 

0.2742 

0.2778 

0.2792 

03<^) 

0.3429 

0.3545 

0.3573 

0.3583 

0.3587 

03'«> 

0.2886 

0.2863 

0.2858 

0.2856 

0.2855 

TABLE  19 


of  Eq.  (5.21)  for  Composite  W,  (90/9') 


0' 

15“ 

30“ 

45“ 

60“ 

75“ 

o/P) 

-272.41 

-59.540 

-19.650 

-6.2877 

-1.0366 

0  •(?) 

1 

-28.239 

-5.8883 

-1.7000 

-0.2970 

0.2543 

c/P) 

-3.9066 

0.7872 

1.6667 

1.9614 

2.0771 

_ 1 

-3.4539 

0.8866 

1.7000 

1.972- 

2.0795 

-.03907 

0.0787 

0.1667 

0.1961 

0.2077 

0.1477 

0.1970 

0.2063 

0.2093 

0.2106 

P3«> 

0.2014 

0.2117 

0.2137 

0.2143 

0.2146 

03'«> 

0.2024 

0.2119 

0.2137 

0.2143 

L. 

0.2146 

TABLE  20 

of  Eq.  (5.21)  for  Composite  T,  (90/0') 


0' 

15“ 

o 

o 

I 

0 

60“ 

75“ 

o 

1 

-228.12 

1 

•61.789 

-21.688 

-7.4310 

-1.7076 

-15.948 

-4.4035 

-1.6200 

-0.6304 

-0.2332 

-2.9311 

0.3289 

1.1149 

1.3943 

1.5065 

.(p) 

2.4089 

1.7733 

1.6200 

1.5656 

1.5436 

-0.5329 

0.0598 

0.2027 

0.2535 

0.2739 

0.2232 

0.2643 

0.2742 

0.2778 

0.2791 

03«) 

0.2695 

0.2812 

0.2840 

0.2850 

0.2854 

«3'"’ 

0.2886 

0.2863 

0.2860 

0.2856 

0.2855 

TABLE  21 


^  (2)  ^  (2)  ^  (3) 

*  °3  ’  ®1 


1  *  “3  ’1  *3  '1  ’  3 

of  Eqs.  (5.21,  5.24,  5.26)  for  Composites, 


(3) 

(0/90) 


(0/90) 

Composite  1 

Composite  2 

Oj^^P^  X  lO^psi 

0.3965 

0.4032 

Oj^'^P^  X  lO^psi 

0 

0 

0^^^^  X  lO^psi 

20.083 

22.113 

a^'^P^  X  lO^psi 

2.1000 

1.5400 

(2) 

"l 

0.1888 

0.2618 

ai-(2) 

0 

0 

0.2497 

0.3533 

03-^2) 

0.2100 

0.2800 

0.1050 

0.0700 

1 

1 

0.0220 

0.0196 
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